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ABSTRACT, CONCLUSIONS AND RECOMMENDATIONS 

The purpose of combining two Quarterly Progress Reports under one 
cover is twofold; first, the breadth of the topic treated here (viz. 
bridging the gap between theory and experiment) is such that the task 
was not completed at the end of the first quarter, and breaking the 
results into two parts would only have caused unnecessary detouring of 
the readers' effort; second, the degree of work intensity was less than 
before since some staff was on summer vacation during the first quarter, 
while a portion of the second quarter was in the period of the no-cost 
extension of this contract, when one of the investigators (Mr. Messalem) 
had left for another appointment. It is of interest that in spite of 
these factors, the progress achieved in the last two quarters seems to 
us of decisive importance. 

A method of evaluation of the transport experiments was worked out 
which is based entirely on conservative fluxes, i.e. fluxes of quantities 
which do not vary across the membrane in the steady state (the "volume 
flux" and "separation flux" which have been used frequently in the past 
by other investigators are not conservative). Using this procedure the 
conductance coefficients for the system 

I * I 

C-103 cation-exchange membrane 0.1 N NaCl 
were calculated. In these calculations two reciprocity relations (one 
relating to salt flow - water flow coupling in osmosis and hyperfiltra- 
tion L = L , the other relating membrane potential to cation transference 
sw ws 

number L = L ) were assumed to be satisfied; the third (relating streaming 
se es' 

American Machine & Foundry Co., Stamford, Connecticut. 
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potential to electroosmosis, L ew = L ) was proven to hold, within 6 per-. . 
cent, at least up to a pressure difference of one atmosphere. In the 
future it should be possible to obtain sufficiently accurate membrane 
potential measurements to determine the range of validity of the reci- 
procity relation L = L . We do not expect to be able to verify the 
reciprocity relation L C11 = L in our apparatus, since the pressures 

b W Wb 

necessary to determine the hyperfiltration coefficient, L , with reasonable 

ws 

★ 

accuracy are higher than our apparatus can support. 

The calculation of the conductance coefficients, L, which characterize 
the transport properties of the membrane, was performed with the experi- 
mental data available from recent and earlier experiments. Not all of 
these data were of the high accuracy which we can now achieve; moreover 
some of the data came from experiments in which the membrane was in con- 
tact with 0.1 N NaCl solutions on both sides, rather than 0.1 N on one 
and 0.05 N on the other; but several importants points emerge from con- 
sideration of the calculated L-values (Table 4) and the methods used to 
calculate these values; viz. (a) The coefficients related to a larger ex- 
tent to water, w, than to salt, s, are larger (L ew » L es ; L ww » l_ ws ; L $w » L ss 
subscript e stands for the electric current). This is probably due to the 
much higher concentration of water than of salt in the ion-exchange membrane, 
for the coefficients L-. increase in general with increasing c-. This property 

I J sJ 

of the set of L-coefficients thus probably reflects the exclusion of salt 
from the ion-exchange membrane ("Donnan effect"). 

★ 

This limitation was known since the principal investigator's early high- 
pressure hyperfiltration experiments performed in steel vessels (1). The 

- advantages- of--a- plasti c -eel 1 -for- the concentration-cl amp-apparatus -seemed 

so large, however, that it was decided at the start of this research to 
forego the determination of the validity range of the reciprocity relation 
L ws = L sw ; the other two reciprocity relationships can be checked in our 
concentration-clamp apparatus. 


2 



(b) The numerical calculations in section A. 2. of the Appendix show the 
large dependence of all L-coefficients on the membrane conductance . In 
fact, membrane conductance emerges as the most decisive transport charac- 
teristic, thus justifying a posteriori the time-honored decision of battery 
technologists to use conductance measurements as primary screening tests 
for separator materials. The implication of this conclusion for future 
experiments is the emphasis on refinements of the technique of membrane 
conductance measurements. A more elaborate and more accurate technique 
than the one used in the past is being developed. It uses the existing 
apparatus described in Progress Report No. 12. Apparatus and procedure 
were modified to yield much more accurate measurements of the distances 
of the electrodes from the membrane, (c) Consideration of the numerical 
values of all terms in the calculation of the conductance parameter, L $s , 
which relates the rate of electrolyte diffusion through the membrane to 
the "electrolyte diffusion force", -Ay s , (determined by the concentration 
gradient of the electrolyte) leads one to the conclusion that osmosis +-*• 
dialysis with shorted electrodes (rather than with unconnected electrodes) 
should be tried, in an attempt to improve the accuracy of this determination. 
Also it would be worthwhile to investigate under what conditions, if any, 
the planned, simultaneous application of several independent forces improves 
the accuracy of the determination of the conductance coefficients, L. 

Finally, the methods developed here for calculating a set 
of L-coefficients from specified transport measurements lead to the unequi- 
vocal characterization of any system of the type: solution | membrane | solution 
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provided only that these specified measurements can be performed with 
sufficient accuracy. The "concentration-clamp" apparatus developed in 
this project represents a large step towards the attainment of this experi- 
mental aim. 

In the future we propose to (A) measure those transport coefficients 
(e.g. electrical conductance) for the system 

0.05 N NaCl | C-103 | 0.1 N NaCl 
which we had previously measured only for the system 

0.1 N NaCl | C-103 | 0.1 N NaCl , 

(B) use the method described in this report for systematic studies of 
other solution-membrane systems and (C) investigate the points mentioned 
under (c) above. 
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L ,L ,L .etc. 
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P 
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T 
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*1 


V i 
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Cl 



A<j> 


o 

effective surface area of membrane, cm 

-3 

concentration, mole cm 

_3 

average between left and right concentrations, mole cm 
membrane thickness, cm 

generalized force (e.g. Ap), joule mole"^ 

_2 

electric current density, amp cm 

-2 -1 

flux of component i, mole cm sec 

-2 -1 

"separation flux" defined in eq. (14), mole cm sec 
volume flux, cm sec - ^ 

conductance coefficient in "Onsager set of transport 
equations" (3)-(5). i, j = w; mole 2 joule -1 cm -2 sec -1 

conductance coefficients in "Duncan set of transport equations" 
(23)- (25) . For units, see Table 1 

conductance coefficients in "Michael i-Kedem set of transport 
equations" (35)-(37), mole 2 joule" 1 cm -2 sec -1 

ic 

pressure, dekabar 

gas constant, 8.31 joule mole ^ (°K) - ^ 
temperature, °K 
time, sec 

transference number of component i, relative to the membrane 

_3 

partial molar volume of component i, mole cm 

3 -1 

= v AgC1 - v Ag [eq. (Al-13)], cm mole 

valency of particle i, equiv mole”^ (positive for cations, 
negative for anions, zero for water) 

5 -1 

Faraday's constant, 0.965 x 10 coul eq 

electric potential difference across membrane, volt 


1 dekabar = 1 joule cm -3 = 10 bar = 9.8692 atm. 

'k'k 

often measured with "Luggin capillaries". 
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A<|>_ potential difference between two Ag/AgCl (anion-reversible) 

electrodes in solutions adjacent to membrane, volt 

k conductivity, ohm"^ cm~^ 

p resistivity, ohm cm 

u.j total potential of component i [eq. (26)], joule mole"^ 


y- chemical potential of component i (includes pressure-volume 

term) , joule mole -1 

Q 

y. concentration-dependent part of chemical potential of 

component i, joule mole -1 


Subscripts' 


+ cation 

anion (Cl~ in this report) 
w water 

s electrolyte (NaCl in this report). The 

equations refer to an electrolyte composed of a cation of 
valency +1 and an anion of valency -1 ("1-1" electrolyte) 


Sign Conventions 


(1) Positive direction is from left to right. (2) Fluxes from left 
to right are counted positive. (3) The operator. A, for finite differences 
refers to the value on the right (double primes) minus the value on the left 
(single primes), as does conventionally the differential operator, d. 

Driving forces are of the general form (-dy/dz). Thus positive values 
of the driving force, (-dy/dz) > 0, lead to positive fluxes. For example. 
Ohm's law is 


i = (1/p) ( -A# ) 

"Driving force" 

and Fick's law 

J i = D i ( ). 

"Dri ving ..force" 
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THEORY 


We consider the isothermal system 

Ag/AgCl | chloride solution, c' membrane | chloride solution, c" | Ag/AgCl 


There are three types of particles which can migrate in the membrane, 

viz. cations (+), Cl” (-) and water (w). Three equations relating the fluxes, 

-2 -1 -1 
J (mole cm sec ) to the generalized forces, F (wattsec mole ) are 

necessary to characterize the transport properties of the system. 


The most fundamental set of transport equations is based on the classical 
concepts of non-equilibrium thermodynamics (2) originally developed by 
Onsager (3) and applied by him to diffusion in liquids (4). Previous exper- 
imental work justifies the description of the fluxes as sums of the products 
of generalized forces and conductance coefficients, L. .: 

' J 


J. = E L. . F. 
1 j TJ J 


( 1 ) 


From a practical viewpoint it is essential to look for transport descrip- 
tions with reasonably constant conductance coefficients (i.e. conductance 
coefficients which are reasonably independent of the forces). Only in this 
case does a single set of L, . coefficients characterize a given membrane- 

1 J 

solution system, i.e. enable us to predict the fluxes from knowledge of 
the forces and/or vice-versa. 


For ion and water transport in membranes under the influence of dif- 
ferences of pressure, concentration and/or electric potential, the generalized 
force acting on particle i is 

F. = v i Ap + Ay!: + z i (2) 


7 



* c -1 

where Ap (dekabar ), Ay.. (wattsec mole ) and Acj> (volt) are the differences 
of pressure, of the "concentration-dependent part of the chemical potential" 

[ i.e. PM Ac. = RTAlna. and the electrical potential across the mem- 

\ 3c i) J 

brane respectively. Hence the original set of transport equations (which 
we shall call the "Onsager set") is 


J + = ~ L +w^w Ap + Ap w> ‘ L ++^ + a P + AU+ + ^A<j>) - L + _(v_Ap + Ay^ - m) 0) 

J _ = " L _w ( V p + ' L -+(VP + AU+ + - L (v_Ap + Ay^ - (4) 

J w = ‘ L WW^W AP + AP W } ‘ L W+^+ Ap + A *4 + " L wJv_Ap + Ay^ - fty) (5) 


Onsager's theoretical work (3) based on statistical mechanics led to 
the conclusion that the conductance coefficients, L-., satisfy the reciprocity 

* J 

condition 

L ij ■ L ji < 6 > 

Onsager's work is of general nature and covers many more phenomena than 
transport in membranes, e.g. heat conduction and chemical reactions. In 
general, the range of validity of linear laws of the type of eq. (1) is 
limited; for instance for most chemical reactions the linear approximation 
is reasonably adequate only very close to chemical equilibrium. If this were 
generally true for transport in membranes also, the linear equations would 
therefore be only of very limited use; for many transport phenomena, however, 
the range of approximate linearity extends far beyond equilibrium. For 
instance, it has been known for many decades that the ion flux in electrolyte 
solutions is proportional to the applied electric voltage over a considerable 

*1 dekabar h 10 bar - 9.87 atm; this pressure unit is used in order to mini- 
mize the use of conversion factors in our system of units (see List of Symbols) 
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voltage range, provided proper electrodes and stirring devices are used 
("Kohlrausch 's law"). Also self-diffusion of particles follows a strictly 
linear relationship between flux and force even in solutions and gases 
which are very far from isotopic equilibrium [Fick's law of diffusion can 
be shown to be a linear law in the sense of eq. (1); see, for instance 
refs. (5) and (6)]. Moreover, many membranes and porous media satisfy 
d'Arcy's law (proportionality between water flux and pressure difference, 

Ap) up to quite high values of Ap. Therefore the range of applicability 
of the linear transport equations (3) - (5) might prove to be appreciable, 
at least with respect to variations of pressure, Ap and electrical potenti al 
differences, A<j> (7). As for Ay C , the situation is more complex. There 
are membranes in which some conductance coefficients are only mildly 
affected by changes in the concentrations of the solutions bracketing the 
membrane. For instance, present knowledge on the nature of ion-exchange 
membranes leads one to believe that changes in the concentrations of 
dilute solutions (up to about 0.5 N NaCl ) should only very moderately 
change the cation conductance coefficient, L ++ in high-capacity cation- 
exchange membranes (such as C-103 used here) .since the "Donnan-effect" 
excludes the anions from the membrane to a large extent. On the other hand, 
a relatively much larger change of the anion transport coefficient, L__, 
with the concentration of the bracketing solutions is expected; this coeffi- 
cient varies very strongly with the concentration, c_, of the anions in the 
membrane (8, eq. 17), and, while of small magnitude, percentage-wise c_ 
varies considerably with the concentrations of the bracketing solutions. 



In other words, the electrical conductivities of these membranes when 
equilibrated with pure water and 0.1 N NaCl vary only by a few percent 
since in both cases the electric current is carried primarily by cations 
whose concentration in the membrane varies relatively little in the two 
cases. On the other hand, anion transport and salt diffusion are controlled 
by the anion concentration distribution in the membrane which varies rela- 
tively much in the two cases, being zero for the membrane equilibrated 
with pure water and small, but finite, for the membrane equilibrated with 0.1 N 
NaCl (7). In non-permselecti ve membranes, both L ++ and L__ vary strongly 
with the concentration. Therefore one can not a priori expect constancy 
of all the conductance coefficients with respect to variations in the 
solution concentrations. (The "concentration-clamp" method was invented 
with this thought in mind; it keeps the solution concentrations constant 

at least while each experiment for the measurement of transport parameters 

* 

is in progress). Since one is interested in an invariant set of transport 
parameters for each system consisting of a given membrane and a preferably 
wide concentration range of solutions of a given electrolyte, attempts 
have been made to break down the conductance coefficients into products 
of concentrations and so-called "friction coefficients" in the hope that 
the latter vary only mildly with the concentration (8,9). Preliminary 
results of this type of analysis by many investigators [summarized in ref. 

(10)] look quite promising; this report does not present this analysis, 
however, since we don't have enough experimental data on our system yet. 

i.e. independent of the concentrations and all forces. 
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The practical application of the Onsager set of transport equations 
[eqs. (3) - (5)] is hampered by the lack of knowledge of the partial 
ionic volumes v + and v (as opposed to the partial molar volumes of 
water and electrolyte which can be unequivocally defined from macro- 
scopic measurements and often found in tables of properties of electro- 
lyte solutions). Moreover, the electrical potential difference across the 

★ 

membrane, A<j>, appears in two of the three forces , whereas the potential 

difference actually measured, Acf>_, is that between two Ag/AgCl electrodes 

placed at some distance from the membrane. Therefore previous investigators 

have performed certain transformations of the Onsager set of equations to 

make it directly amenable to evaluation of the experimental results. We shall 

critically discuss two transformations, one by Duncan (11) who initiated 

this project, and one found in the text of Katchalsky and Curran (2). 

We shall then justify the use of a modification of the latter transformation 

** 

for the evaluation of our measurements and demonstrate this evaluation. 

The Duncan Transformation (11) 

The purpose of this transformation is the use of driving forces which 
can be readily measured instead of the generalized forces in the original 
Onsager equations (3) - (5) which involve ionic volumes and the electrical 
potential difference across the membrane. The "price" paid for this simpli- 
fication of the forces is the introduction of two non-conservative fluxes 
(i.e. fluxes which are not uniform across the membrane even in the steady 
state). The impact of this problem will be discussed after presenting 

•k 

— It-i s-not-easy to measure -A<j>-accurateTy-when a- current-passes -through 

the membrane. 

This transformation is referred to as the "Michaeli-Kedem" ("M-K") trans- 
formation in this report. 


11 


the transformation in elementary manner. In the following derivation 

consider a membrane bracketed by two solutions of the same electrolyte 

With a concentration difference, Ac, so small that in the steady state 

the difference between the volume flux, Jy , leaving one solution to 

enter the membrane at one face, and, after passage through the membrane, 

* 

entering the other solution, is negligible . (A similar consideration 
is to hold with respect to the "difference flux", J Q , to be defined 
later). With this reservation, it is possible to define a volume flux, 
Jy (cm cm sec" = cm sec ) without specifying the location at which 
it is considered: 


J v 
w w 


J + v + + J. 


(7) 


The i/.'s are the partial molar volumes of the components; one takes 
the average of their values in the two solutions bracketing the mem- 
brane. Substituting for the fluxes J + , J_ and J w from equations (3) - (5) 
respectively, we obtain: 



L VV 


-V-ww'V + L + v/+ + L V-> - ^ L w + \ + L ++*+ + L - + vJ + ^_ C t L w_v w + L +-*+ + 


-S^fl-w+Vw + L ++ v + + L_ + v_ - L w _v w - L + _v + - L__v_) 


( 8 ) 


-VE 


At first sight, it might seem as if in the steady state the volume loss of one 
solution were numerically equaT to the- volume-gain of- the others -Inasmuch as -the-- 
partial molar volumes of water and electrolyte differ somewhat with the solution 
concentration, volume conservation does not prevail, however, i.e. the volume loss 
of one solution does not quite equal the gain of the other especially if Ac is 
appreciable. 
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J v = L vv^ 

- Ay 


Now transform the center term of eq. (8) which contains the chemical 

potential differences, Ay C , by using the Gibbs-Duhem relation [ref. 2, 

eq. (9-4) to substitute for A y^ in terms of the chemical potential of 

the salt, Ay c : 
s 


1 Ay f f 


(9) 


s / T c w 

(valid for small concentration differences between the two solutions). 

E s and £ w are the averages of the salt and water concentrations in the 
two solutions respectively. Also substitute for the potential difference 
across the membrane, A<}>, in terms of the measured potential difference 
between the Ag/AgCl electrodes, A<j>_ , from eq. (AT -17): 


- W = - ^A(j>_ - Ay_ 


( 10 ) 


This transformation yields 


e 


•AP) + Au.(L v ^ * L_5 + i + L ! h. ) - iP + C (L w+ v w + L ++ v + ♦ L. + 5_) 


s' ww w - +w 

c c 

w w 


w 


+ L + . ; + + L„v.) + L VE ( -?AO - AP-O-w^w + >■++*+ + L_ + v_ - L w .V L + -V 

( 11 ) 

The chemical potential of the salt, y^, is the sum of the chemical 


potentials of its constituent ions. Therefore 

Ay g = Ay+ + Ay^ 


( 12 ) 
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= L vv^ 


Hence eq. (11) reduces to 

■ i(,) + f L «4 + L ++ '' + + L - + ; - - T ( l „w 5 w + L ^ + + L - w i '->] <-^ s C) + L VE (-5 t 'M.) 

(13) 


w 




■VD 


The second flux used in the Duncan transformation is the "separation 

O -1 

flow", (mole cm sec ) which is proportional to the velocity, u s , of 
the salt relative to that of the water, u^ (ref. 2, p. 99). 


11m °D = J $ - c- J w = c s (u s ■ u w> 
AC— >0 W 


(14) 


In the steady state, the fluxes J g and J w are uniform, but since the 


concentrations vary across the membrane, Jp, u g and u w are not uniform 
("conservative"). Again, for small concentration differences this non- 
uniformity is neglected and average values of the parameters are used: 
c_ 


J D =J s-f J w 
c w 


(15) 


In the experimental system considered here the salt flux, J s , is equal 

in magnitude to the flux of sodium ions, J + (both being measured in mole 
-2 -1 

cm sec ), because to maintain electroneutrality, the Ag/AgCl electrodes 

supply and subtract an equivalent amount of Cl for the amount of Na + 
entering or leaving a compartment respectively: 

J s = J t (16) 

Substituting J + for 0 S in eq. (15), and then substituting the expressions 
for J + and J w from the original flux equations (3) and (4) respectively, we 
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obtai n 


J D = J + - y J w 
c w 


-A P [( L +w v w + L ++ v + + L + _v.) - Lww v w i - L w+ v + ^ - L w _v_ ^ ] 




~^DV~ 


w 


w 


'w 


‘ Ay w^ L +w “ - L ww^ " Ap +^ L ++ " - L w+^ 


- AP. C (L + . - I 1 L w _) 


w 


w 


'w 


c c c 

Si i _i_Si \ a.C/i s 


•^A^_(L ++ - ^ L w+ - L + _ + ^ L w .) - Ay‘(L ++ - ^ L w+ - L + _ + ^ L 

V c w ^ c w J c w c w 

L DE (17) 


w 


Combining the terms containing Ay+ with those containing Ay c by use 
of eq. (12) we obtain 


J D *~D v^" A P^ + ("^s) ^ L ++ " - *~w+^ + (~ Ay J^ L +w ” - ^ww^ + W' 3 ^ 


w 


l 1 


w 


Substituting for Ay^ in terms of the chemical potential of the salt, 

c * 

Ay s> by use of the Gibbs-Duhem equation (9), we obtain 

c c c 2 

J 0 - L DV (-AP) + (-Ay') (L ++ - ^ L w+ - L + „ ^ + * L( J + L pE < - 3 Ta<|,_ ) (19) 


w 


c w c w 


-DD 


_? 

Finally, the current density, i (amp cm ) divided by Faraday's constant 
(coul eq ^ ) is the third flux in Duncan's set of flux equations. In the 
steady state, this flux is conservative, because of electroneutrality. By 
using the expressions for J + and J_ from the original flux equations (3) and 
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-(4) we obtain 

i/^' = J + - J - = (“ A P) ( L +v /w " L Vw) + (L ++' ? + ' L -+*+) + ( L +-V. - L -_vJ] 

V. 


"^w^w - L -w> - A ^ L ++ - L - + ) - A P-< L + - - L ~> 


■^A(j.(L ++ - L_ + - L + _ + L__ ) 


( 20 ) 


’EE 


Substituting for -Ay^ from the Gibbs-Duhem equation (9) and for 


■ A<}> from (10) we obtain 


c c 

c, , s , , s 


1/y - L EV (-Ap) - A|£(-L ■+ L “ ) Au^(L +t - L_ + ) - Ap^(L + , - L„) 

c, 


w 


w 


+ L EE (-&A<p_) - Ay_ (L ++ - L_ + - L + _ + L__) 


( 21 ) 


p P 

Combining the terms containing Ay + and Ay_ by means of eq. (12), we 
obtain 

i/y - L EV (-Ap) + (L ++ - L + + L_ w L +w i ) (-Ah') + L EE (-? A*. ) (22) 

c w c w 



It is seen that if the reciprocity relations [eq. (6)] are satisfied 
in the original flux equations (3) - (5), i.e. if L + _ = L_ + , L +w = L w+ , 

L = L , then there is also reciprocity in the transformed equations, i.e. 

L VD = L DV’ L VE = L EV’ L DE = L ED* We cor,clude therefore that in the transformed 
equations, the fluxes Jy, J D and (i/^) are properly conjugated to the driving 
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£ 

forces -Ap, -Ay s and - 3^A4>_ respectively. In summary the three transformed 
flux equations (13), (17) and (22) are rewritten here: 


+ Lyp(-Ap s ) + Ly^(-^A(J)_) 


(23) ■) 


(24) 


I "Duncan set 
^ of flux 
f equations" 


i/^ = L EV (-Ap) + L ED (-Ay^) + L EE (-9 r A4>_) (25) 

J 


It should be noted that the equations given by Duncan in ref. (11b) 
are arranged in somewhat different form; specifically in the last equation 
- yA<f>_ appears on the left and i on the right; also the terms containing 
the driving force -^Acf) in equations (23) and (24) have been replaced by 
terms containing the conjugated flux, i/^ , instead. Since both the driving 
forces and the fluxes in eqs. (23) - (25) are almost identical with those in 
ref. (11b), however, we feel justified in applying the term "Duncan set of 
flux equations" to eqs. (23) - (25). 

Table 1 lists the conductance coefficients, L, in the Duncan set of 
flux equations in terms of the original coefficients in the Onsager set, 
eqs. (3) - (5). 

The relative simplicity of equations (23) - (25) derives in part from 
the fact that they contain the measurable potential difference, A4>_, between 
the Ag/AgCl electrodes rather than the potential difference, AcJ>, across the 
membranes which can only be calculated if assumptions are made about the 
difference of individual electrode potentials or junction potentials (14). 
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Relations Between Conductance Coefficients in Onsager-Type Set of Flux 
Equations (3) - (5)* and Transformed Set (B. C. Duncan), Eqs. (23) - (25) 
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On the other hand, the use of the non-uniform fluxes Jy and J D limits the 
direct utilization of these equations to systems with very small concentra- 
tion differences between the two solutions bracketing the membrane. 

The Michael i-Kedem Transformation (2, p. 151, 15) 

This transformation leads to three fluxes, J , J w and i/^, which are 
all uniform in the steady state. Moreover the three conjugated forces, 
|-Ay s |,(-Ay w | and Ay_ = -3»A<j>_ + (v^-j - v Ag ) Ap can be readily measured 
and/or found in thermodynamic tables. For this reason we have now come to 
the conclusion that this set of equations is the most practical for the 
evaluation of our past and present transport experiments, and that the 
L-conductance coefficients in this matrix are the most practical for char- 
acterizing the transport properties of the membrane. [It should be noted 
that we consider this type of characterization as an important intermediate 
step, but plan to eventually use friction coefficients (8 - 10) for mem- 
brane characterization]. 

The starting point for this transformation is the original (Onsager) 
set of flux equations (3) - (5) which is rewritten, using the symbol Ay i 
for the difference between the total potentials of component i in the right 
and left solution respectively [see also eq. (Al-1) in the appendix]: 

o 

Ay.. 5 Ay i + = V^Ap + Ay^ + Z^Acp (26) 
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Since the net charge of water, z w , is zero Ay w = Au w . 

With this notation, the set of original Onsager flux equations is 

J + = L ++ (-Ay + ) + L + _(-Ay_) + L +w (-Ay w ) 


j_ = L_ + (-Ay + ) + L C-Avi_) + L (-Au ) 


J = L , (-Ay.) + L (-Ay ) + L (-Ay ) 
w • w+ ' M + ' w- v ' ww ' K w ' 


(27) 

(28) 
(29) 


The total potential of the salt component, Ay s , is equal to the sum 
of the total potentials of the constituent ions. Also, since salt carries 


no net charge, z = 0, and hence Ay = Ay c 


Ay s = Ay + + Ay_ 


We substitute Ay $ - Ay_ for Ay + in eqs. (27 ) - ( 29 ) and then 
-3>A<f> + Vq-j A p for Ay_ from eq. (Al-14) in the appendix. The result 
is the following set of equations: 


(30) 


J + -J s - L ++ (-Ay s ) + (L ++ - L + J (-^AcD_ + V^Ap) + kJ-Ay*) 


(31) 


J_ = L_ + (-Ay s ) + (L_ + - L ) (-# A4>_ + v^Ap) + L_ w (-Ay w ) 


J w = L w + (- A ^s ) + (L ^ - L w-> + V C1 A P) + L ww (_Ay w ) 


’w+ w 


cr 


ww' W' 


(32) 

(33) 


Finally, we calculate i/^ from eqs. (31) and (32): 

=■ j + - o. - (L ++ - L_ + )(-au s ) + (L ++ - L t . - L_ + + + » c ,Ap) + < L « ‘ L -w )( - A|J w ) 

(34) 
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Summarizing, the transformed set of flux equations (31), (33) and 
(34) is: 


J s = L ss ( ‘ Ap s ) + L sw ( ‘%> + L se (-S' 4 *- + v cl Ap) ' 

J w = L ws ( ' Aw s ' + W’V + + v Cl Ap > / 

V? = L es (-Ap s ) + L ew (-A Uw ) + L ee (-3rA*. ♦ v cl Ap) , 


"Modified M-K 
set of flux 
equations" 


(35) 

(36) 

(37) 


where the meaning of the new L conductance coefficients in terms of those 
of the original set, eqs. (3)-(5) is shown in Table 2. Note that when 
reciprocity [eq. ( 6 )] is satisfied in the original set, reciprocity also 
prevails in the transformed one. We conclude that in eqs. (35)-(37) fluxes 
and forces are properly conjugated. 


Experimental Results 

Since the "concentration-clamp" apparatus and the operating procedure 
have been described in detail in previous reports and in a recent publication 
from this project (16), no details about these aspects of the work are 
presented here. The cell for the measurement of the membrane conductivity 
and the measurement procedure have been described in Quarterly Report No. 12, 
Section A. III. 4 (p. 15-17) and therefore they are not repeated here either. 
The results used for the calculation of the L-coefficients are listed in 
Table 3 . In addition, we shall use the value of 5.7 ± 0.8 ohm cm^ for the 

p 

resistance of one cm of active (uncovered) membrane area. This value 
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TABLE 2 


Relations Between Conductance Coefficients in Qnsager-Type Flux Equations 
( 3 ) - ( 5 ) and Transformed Set (I. Michaeli and 0. Kedem) Eqs. (35)-(37) 


ss 

= 

L ++ 


ww 

= 

L ww 


ee 

= 

L ++ - 

L + - 

sw 

= 

L +w 


ws 

= 

L w+ 


se 

= 

L ++ - 

L + - 

es 

= 

L ++ - 

L - + 

we 

= 

L w+ " 

L w- 

ew 

= 

L +w 

L -w 


2 -1 -2 -1 

All coefficients have the units mole joule cm sec 



TABLE 3 


Experimental Results Used to Calculate M-K Conductance Coefficients, L 


Experiment No. 

163 

159 

150 

145 


Permeation under 

Electromigration 

Dialysis-*-^- 

Type 

pressure 

- electroosmosis 

osmosis 

Pressure difference, Ap, atm 

-1 

-0.8 

0 

0 

★ 

Salt concentration on right , 

c" , mole cm -3 
s 

10 -4 

1 

O 

10“ 4 

1 .001 6 x 10" 4 

★ 

Salt concentration on left , 
c^, mole cm -3 

'3' 

1 

O 

10' 4 

1 

O 

1 

O 

i — 

X 

LO 

CTi 

o 

Concentration difference, Ac , 
mole cm -3 

0 

0 

0 

0. 506 o x 10" 4 

?** 





Current density, i, ma cm~ z 

0 

0 

-0.971 

0 

*** 6 





J v cm sec -1 x 10 

Salt flux, J ,mole cm -2 sec -1 
x 10 9 s 

0.374 

0.304 

-1.805 

0.284 



-8.69 

-0.0713 

Water flux, J. mole cm -2 sec -1 





x 10 9 w 



-99.79 

15.3 

Potential difference, A<f> , volt 





x 10 6 





Cation transport number, t + 



0.864 


Water transport number, t 

w 

Streaming potential, A<f> , volt 



9.63 


x 10 3 

0.174 

0.146 




Buret side, right; column side, left; see Figure 1. c determined by chloride titration. 
Positive electrode on right side. Active membrane area 8.13 cm . 

"jc'fc’Jc 

Calculated from data on the right (buret) side only. 
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'FIGURE l. SCHEMATIC REPRESENTATION OF THE CONCENTRATION FEEDBACK MECHANISM 




corresponds to recent results obtained in the conductance cell described 
i,n Progress Report No. 12. [The active area of the membrane (not covered 
by perforated supports) was 20.2 cm and the membrane thickness 0.017 cm]. 

A recent critical reexamination of the raw data from the measurements, 
described in Progress Report No. 13, and observation of a progressive change 
in the color of the active' areas of the membrane made it seem advisable to 
remeasure the resistance. In the process of remeasurement, the procedure 
was refined and additional measurements by this new procedure are in pro- 
gress. (It should be noted that exact measurements of the resistance of 
thin, conductive membranes are by no means a simple matter.) Results of 
the newest measurements evaluated, which will hopefully be more accurate, 
should be available soon. 

Method for Calculation of the L-Conductance Coefficients for the M-K 
Transport Equations 

In the strict sense the Michaeli-Kedem ("M-K") transport equations, 

(35 ) - (37 ) apply to the total system between the electrodes, i.e. the mem- 
brane plus the two adjacent solutions. Since we are primarily interested 
in the transport coefficients of the membrane proper, we inquire first 
which of the fluxes and forces would be appreciably affected if we progressively 
reduced the thickness of the solution layers in series with the membrane. 

Since the solutions are wel 1 -stirred, concentration gradients in them 
are negligible. Thus the generalized forces -Ay c and -Au which represent 
the chemical potential differences between the solutions near the electrodes 
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are identical with the corresponding differences across the membranes 
only. As for the electrical potential difference, A<j>_, it will depend 
on the magnitude of the current density, i. When the latter is zero, 

A<f>_ will be independent of the position of the electrodes, and hence 
the experimental A<j>_ measured at a finite distance, can be substituted 
for the corresponding potential difference of electrodes in very close 
vicinity of the membranes. When the current density is finite, however, 
and kept constant, moving the electrodes towards the membrane will 
decrease the absolute value of Acj>_, since the ip drop through the solu- 
tions is decreased. Therefore, if we want to calculate transport coeffi- 
cients for the membrane only, without contribution of the solution layers, 
we have to make allowance for this ip drop, i.e. instead of using the 
measured values of A<j>_, we use the value corrected for infinitesimally 
small distance from the membrane: 

(Ad> ) . = (A<j) ) . + i ( p ' z ' + p"z") (38) 

where the p's are the specific resistances of the solutions and the z's 
the distances between the membrane surfaces and the electrode surfaces 
in the two solutions respectively. 

This is exactly what our method of measuring the membrane conductance 

I 

* 

in the cell with the movable electrodes makes possible. Moreover, when 
the current is negligible (e.g. in measurements of streaming potentials 
or membrane potential s) ,it is seen from eq. (38) that the measured potential 
differences are independent of the position of the electrodes. 

_ 

The cell is described in Quarterly- Report- No. 12 
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In a thought experiment one can reduce all the L-coefficients 
to a system in which the electrodes are so close to the membrane sur- 
face that the mass and ion-transfer resistances of the remaining thin 
solution layers are entirely negligible. The L-coefficients calculated 
for this system may be considered as the transport coefficients for the 
membrane. 

The methods of calculation of the conductance coefficients from 
the experiments described in Table 3 are presented in this section. 

The numerical calculations are in Section A. 2. of the Appendix. 


When the solutions are of the same concentration and at the same 
pressure, the algebraic sum of the two Ag/AgCl electrode potentials 
(i.e. of the potential drops between each silver metal and the solution 
in contact with it) is zero. Therefore in this case A<J>_ = A<j>. Since 
the conductivity of the membrane is defined as 

k = (-i/A<j>)d (39)* 

if follows from eq. (37) that 

(L ee ) A M =0 = = -i/(A^ 2 ) = k/ ( 3^ 2 d) = (3< 2 pd) _1 (40) 

In the determination of the membrane resistance (per unit active 
area), pd, the solution resistance is eliminated by extrapolation to 
zero solution thickness. Hence ( L ee ) Ay=0 is a transport characterization 

*The negative sign derives from the convention that for positive potential 
difference (<f>" > 4> ' ) , the current. direction (always taken as the flow 
direction of positive carriers) is negative.- _ 
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constant of the membrane alone. In this measurement, platinized platinum 
electrodes carrying alternating current were used, but for Ay=0 this makes 
no difference, since in this case, the electric potential difference be- 
tween these electrodes is the same as A<j>_. 



Dividing eq. (35) by (37) for Ay=0, we obtain 

L se /L ee " tV'^V 0 = V < 41 > 

Hence 

^ L se^Ay=0 = M L ee^Ay=0 

(^s ) Ay=0 experimentally determined in an electromigration - electro- 
osmosis experiment (e.g. No. 150, Table 1). The calculation of (*- ee )^ _q 

from a membrane resistance measurement has already been described [eq. (40)] 
Hence L se can be calculated from eq. (42). 



Dividing eq. (36) by (37) for Ay=0, we obtain 
L we /L ee = C*V^^O^Ay=0 = *w 

Hence 

^ L we^Ay=0 ~ t w^ L ee^Ay=0 


(43) 


(44) 


Thus (L ). n 
' we , Ay=0 

osmosis experiment, 
of L se . 


can be determined from an electromigration - electro- 
simul taneously _ with and-simi -lar 1 to the determination 
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'es 


and L, 


ew 


In eq. (37) we substitute for Ay w from the Gibbs-Duhem equation (9), 
and setting 

Ay w = Ay w + V p ( 45 ) 

[from eq. (26), with z = 0], we obtain the following equation for the 

★ 

electric current density 

1 = 3r[a es 5 s + L ew 5 w ) + < L es ' L e „ ^ > (‘ A ^> ' L ee (y A *- + A P V C1 )] 

c w 

(46) 


This equation can be used to evaluate measurements of the streaming 
potential, which are usually performed with an electrometer or other 
device which draws extremely little electric current. Setting i = 0 
we obtain from eq. (46) 


AJ . _ L es " L ew^s^Sp A c , 1 

Ap s + § 




ee 


L v + L v 
es s ew w 


'ee 


- v 


Cl 


Ap 


(47) 


Thus a plot of the potential difference between the two Ag/AgCl 
electrodes in the two half-cells, -A<j>_ = <}>_[ - <j>" versus the pressure 
difference, Ap, should be a straight line. In general, the two unknowns, 
L , and L „ , can be determined by solution of the simultaneous equations 
for the slope ("streaming potential differential, SPD") and the ordinate 
intercept ("membrane potential, MP"): 

★ 

Since the Gibbs-Duhem relation is not exactly satisfied for finite con- 
centration differences (see footnote p. 49)_,we_ plan, to use Ay and Ay 
separately in the future without relying on this relation for 5 the 

calculation of L and L . 

es ew 
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8(-A<f>_) 

3(Ap) 


SPD 


1 


/ L v + L v 
i es s ew w 


'ee 


- v 


Cl 


to- Vo 2 MP * 


'es 


- L 


'ew 




< E s /E w> . c 
— 4U S 


ee 


(48) 


(49) 


At the present stage of this research, we have reliable streaming 
potential lines only for a system without concentration gradient (i.e. 

r 

Ay s = 0). Such lines are shown in Figs. 2 and 3. Inasmuch as the 

membrane potential is zero these data can yield only L v + L v , 

es s ew w 

rather than l_ es and L gw separately. Therefore, for the purpose of this 

interim report, we have assumed (rather than proven) that the following 

* 

reciprocity is satisfied, as the theory demands : 



(50) 


Since l_ se could be calculated from the measured transport number 
of the cation [eq. (42)], L gw could then be calculated by rearranging 


eq. (48) 


! 3(-A$_) 

(L ewWo = ; w [L ee (3r ^ _ 


+ V C1 ) 


‘-se'y 


(51) 


★ 

This amounts essentially to assuming that the cation transport number 
determined by an electromigration - electroosmosis experiment agrees 
with that calculated from a suitable measurement of the membrane poten- 
tial (17). This assumption has been made long ago for many membranes 
[see, for instance, (18)], although the relatively small contribution 
of the water transport to the membrane potential has not always been 
duly taken into account. 
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0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

PRESSURE DIFFERENCE, Ap, atm 


FIGURE 2. Streaming Potential Across C-103 Membrane (Experiment 157) 

Solutions 0.1 M NaCl . Short vertical and slanted lines connecting 
points referring to same pressure setting indicate hysteresis effect 
after removal and re-application of pressure. All points corrected 
for small asymetry potentials between Ag/AgCl electrodes separated 
by the membrane, observed at the time of pressure release. 
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ELECTRICAL POTENTIAL DIFFERENCE 
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(52) 


< L ewVo * 0/; w )[?L ee ( s PD ♦ ) - L se 5 s ] 


In this manner we tested the reciprocity relation 


( L ewWo = ^ L we^Ay=0 ( 53 ^ 

and found it to hold within six percent over the limited pressure 
range investigated (up to Ap = -1 atm). The numerical calculations 
are presented in Appendix 2. 

In the future we plan to measure streaming potential curves (in- 
cluding measurements at Ap=0, i.e. membrane potential measurements) for sys 
terns with a concentration gradient. The assumption L se = L can then 
be tested experimentally. 


L pp , L, . L ,, and L 
ss ww sw ws 


These remaining four coefficients can in principle be calculated from 
the following four ratios of flux to driving force: 


(^w' /ilJ s)i=0,Ap=0 

Osmotic water flux 

(osmosis -*-*■ dialysis experiment) 

(J s /Ap s )i = o^p_o 

Dialytic salt flux 

(osmosis dialysis experiment) 

(J */ Ap) i-0 

Water flux 

(hyperfiltration experiment) 

(J s /Ap) 1=0 

Salt flux 

(hyperfiltration experiment) 


The first three of these measurements can be performed in the present 
"concentration-clamp" apparatus, but the accurate determination of the 
salt flux due to a pressure difference is not possible in this apparatus 
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which for a number of good reasons was built from a translucent poly- 
carbonate plastic. The latter is not mechanically stable at pressures 
high enough to measure the salt-filtering effect in hyperfiltration 
("reverse osmosis") with a reasonable degree of accuracy . Therefore 
we assume, rather than prove, that l_ sw = L ; as a result we have only 
three unknowns to be determined from the first three measurements. 

The relevant equations are derived in the following. 

For i=0, we rearrange eq. (37) as follows: 


-(•^A4> - v r ,Ap) = 


'es 


qI'-'h/ - r — Ay 
ee 


l ew a 

s + cr ^ 

ee 


w 


(54) 


Applying this equation to the calculation of the sal t flux in a 
dialysis -*-*• osmosis experiment (where i=0) from eq. (35) in which we 
substitute -(v^Ap + Ay?) = -Ay^ for the first two driving forces, 
we obtain 


*Vi=0 


L ss ( -VP 




+ L (-v Ap 
sw v w y 


*£> * [L es ( Ve + 

ee 


Ay C ) + L (v Ap + Ay C )l 
H s' ew v w K H w ;j 


(55) 


Since the dialysis +■*■ osmosis experiment is performed at uniform 
pressure, this reduces to 


Ap=0 



L se L es 



+ 



es^ew^ 

< 

*"ee i 

Ay s c 


(56) 


The principal investigator has at his disposal a complete tubular hyper- 

- filtration apparatus in which this_ salt rejection effect is being measured 

simultaneous with the hyperfiltration fate and with streaming potentials - 

across modified cellulose acetate membranes at pressures up to 100 atm (19). 

This apparatus is in use for another project, however, and would have to be 
modified to provide the data necessary here, and this could not be done in the past. 
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The determination of l ars and L has already been described-- 
see eqs. (40) and (52) respectively. L gs could not be determined by 
experiments, but has been assumed equal to L se [eq. (42)] in view of 
the reciprocity relation (50). The salt flux, ( j s )-| = q ^-q was 
determined in the dialysis experiment. Hence the only two unknowns 
in eq. (56) are L s$ and L $w . 

In the same manner we obtain from eq. (36) for the water flux: 


i=0 


L ws ( -VP 


Ay ) + L (-V Ap 
M s' ww' w F 


Ay ) + (L /L )(L Ay + L v Ap + L Ay C + L v Ap) 
H w y ' we' ee' y es M s es s K ew M w ew w 


C57) 


Since the pressure is uniform, this reduces to 


J w ' 

_ , L we L ew , i 

(. L we L es ^ 

A c 
Ads 

r Ap wl 

ww ‘ . + | 

i=0 L ee 

l ws Li 

\ ee 1 

A C 

Ap w 


(58) 


The unknowns in this equation are and l ws . 


To evaluate the hyperfiltration experiment, we divide eq. (57) by -Ap: 


I M 
i-*”, 


= L v + 
ws s 


i=0 


. c . c 

Ay s - 
— - + L v + — — 
i ww w 
Ap Ap 


(L /L ) (L Ay s + L 
v we y ee /v es 

Ap 


v + L 
es s ew 


Ay, 

Ap 


— + L v ) 
ew w y 

(59) 


The water flux, J l( , is measured as a function of the pressure difference, 

W 

-Ap. The unknowns in this equation are L c and L . 

W5 WW 

Thus it is possible to determine the three unknowns L^, L ws = L sw and 
L $$ by solution of the three equations (56), (58) and (59), provided all 

*Note that the method of calculation for L wg has already Been described [eq-. (44)].- 
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V.' 


experiments were performed under the same concentration gradient. 

Our hyperfiltration experiments were performed at conditions of 
uniform concentration (i .e. y£ = 0 = no osmosis) by "clamping" 
the concentrations of the two solutions adjacent to the membrane at 
equal concentrations of about 0.1 M KC1 . For this case eq. (59) 
reduces to the simple form 

[ V< A P^i=0 ^ws^Ap=0 ^s + ^ww^Ay=0 ^w ” ^we^ee^Ay=0 ^es^Ay=0 + ^ew^Ay=0 

Ap=0 

(60) 

In principle it is not permissible to solve the triplet of 
equations (56), (57) and (60) for L^, L w$ = L sw , and L $s , because the 
first two equations refer to a dialysis experiment which, by definition 
requires a concentration gradient, while the last equation refers to 
an experiment at uniform concentration. In the future, we plan to main- 
tain the finite concentration difference in all experiments (as we have 
indeed occasionally done in the past), but for the purpose of this report 
we shall utilize the results of the latest measurements, as presented in 
Table I, in which the water flux under pressure was measured under con- 
ditions of Ay=0. Since L is similar to a d'Arcy permeability constant, 

ww 

and since it is known that these "constants" vary but little with concen- 
tration (provided the concentration is below about 1 M) , we do not think 
that the value of L would have been appreciably different if the concen- 

WW 

trations of the solutions in the hyperfiltration experiment had been identi- 

* 

cal ynth those in the dialysis «->■ osmosis experiment . In fact it 

*The water flux, as opposed to the coefficient would change very appre- 
ciably, however, since, in addition to the pressure force, the osmotic 
force would drive the water through the membrane. 
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is usually tacitly assumed that the conductance coefficient L, is 
independent of the concentrations of the solutions bracketing the 
membrane (1). As for the coefficient L ws = L sw , some dependence on 
the solution concentrations is expected (20), but this dependence is 
neglected in this report. Consequently, the three coefficients L , 

L 11C . = L„,, and are determined by solution of the simultaneous 
equations (56), (58) and (60). The calculations are presented in 
part A. 2. of the Appendix. The results of the calculations of all 
conductance coefficients in the M-K matrix are presented in Table 4. 

Conclusions 

While it is not easy to draw immediate conclusions on the transport 

* 

mechanism from the phenomenological conductance coefficients, L , several 
important points emerge from consideration of the L-values in Table 4 and 
the methods used to calculate these values: 

First, the coefficients related to a larger extent to the solvent 

than to the salt are larger (L , » L _ ; L » L 1C ; L C11 » L cc ). This 

3 ew es ww ws sw ss 

is probably due to the much higher concentration of water than of salt 
in the ion-exchange membrane, for the coefficients L.. increase in general . 

* J 

with increasing c, [while the flux equations in ref. (8) are not identical 

J 

with equations (35)- (37 ) used here, some insight into the physical meaning 

i 

of the L-coefficients, and in particular the above conclusion about their 
★ 

On the other hand, friction coefficients , which can be calculated from the 
complete group of L-coefficients are useful for this purpose (8-10). We 
plan to perform such calculations in the future. We expect that it will 
~ be~useful to invert- the flux-force equations- first, -i .e. present -the -forces _ 
as sums of terms linear with respect to the fluxes, and calculate the friction 
coefficients from the resulting resistance coefficients, j (2 ) . 
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TABLE 4 


Phenomenological Coefficients for System 0.05 N NaCl C-103* 0.1 N NaCl 


These conductivity coefficients, L, refer to the flux equations (35)- 

-12-2 -i** 

(37). All coefficients have the units joule mole cm sec 


L 

L 

L 

L 

L 


ee 

se 

we 

ew 

ws 



1.89 
1.64 
1.88 
1.77 
1 .73 
1.42 
1.31 


x 10 
x 10 
x 10 
x 10 
x 10 
x 10 
x 10 


]*** 

- 10 *** 

- 10 *** 

-10 

-11 

-8 


L„ assumed equal to L 
es M se 

> reciprocal pair 
L gw assumed equal to L w$ 


★ 

Cation-exchange membrane made by American Machine and Foundry Company, 
Stamford, Connecticut. 

** 3 

1 joule = 1 wattsec = 1 dekabar cm . 

'Jc'ic'ic 

Determined from measurements on system without concentration gradient, 
0.1 N NaCl | C-103 | 0.1 N NaCl. 
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variation with c. (borne out by the results) may be obtained from that 

J 

reference; see eqs. (17) or (32) of reference (8)]. This property of 
the set of L-coefficients thus probably reflects the exclusion of salt 
from the ion-exchange membrane ("Donnan effect"). 

Second, one of the three reciprocity relations was tested by the 
results of our experiments, and was confirmed within about six percent 
for the conditions used in these measurements. 

Third, the numerical calculations in section A. 2. of the Appendix 
show the large dependence of all L-coefficients on the membrane conductance 
(Figure 4). In fact, membrane conductance emerges as the most decisive 
transport characteristic, thus justifying a posteriori the time-honored 
decision of battery technologists to use conductance measurements as pri- 
mary screening tests for separator materials. The implication of this 
conclusion for future experiments is the emphasis on refinements of the 
technique of membrane conductance measurements. A more elaborate and more 
accurate technique than the one used in the past is being developed. It 
uses the existing apparatus described in Progress Report No. 12. Apparatus 
and procedure were modified to yield much more accurate measurements of 
the distances of the electrodes from the membrane. In case future measure- 
ments yield a somewhat different conductance value, al 1 L-coefficients in 
Table 4 will be different, but their ratios will not change much. 

Fourth, consideration of the numerical values of all terms in the 
calculation of the conductance parameter, L $s , which relates the rate of 
electrolyte diffusion through the membrane to the "electrolyte diffusion 
force"-, _-.Ay , ^determined by the .concentration gradient of the electrolyte) 
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DISTANCE, cm 

FIGURE 4. C-103 Membrane. Evaluation of Resistance (25°C). ' 

2 

Solution 0.1 N NaCl . Active area = 20.25 cm . Abscissa repre- 
sents distance from point of closest approach of electrodes to 
membrane (upper line) or electrode distance (lower line). Re- 
sistance of membrane is determined from ordinate distance of 
the two parallel lines. For lower line ordinate intercept rep- 
resents circuit impedance. . 
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leads one to the conclusion that osmosis ■*-*- dialysis with shorted electrodes 
(rather than with unconnected electrodes) should be tried, in an attempt 
to improve the accuracy of this determination. Also it would be worth- 
while to investigate under what conditions, if any, the planned, simul- 
taneous application of several independent forces improves the accuracy 
of the determination of the conductance coefficients, L. 

Finally, the methods developed here for calculating a set 
of L-coefficients from specified transport measurements lead to the unequi- 
vocal characterization of any system of the type: solution | membrane ) solution, 

provided only that these specified measurements can be performed with suf- 
ficient accuracy. The "concentration-clamp" apparatus developed in this 
project represents a large step towards the attainment of this experimental 
aim; the following next steps are desirable: (a) further refinements of 

the apparatus and procedure in order to obtain this desired accuracy in a 
few types of transport measurements in which we still need these improve- 
ments, and (b) the systematic collection of transport data. 
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APPENDIX 


A. 1 . Relation of electrical potential difference, A4>, across the membrane, 
to potential difference, between two Ag/AgCI electrodes. 


Consider the isothermal system 


Ag/AgCI 


NaCl solution, c' , p 1 


Membrane 


NaCl solution, c", p" 


Ag/AgCI 


The condition for equilibrium between each solution and the Ag/AgCI 

* 

electrode in contact with it is that the sum of the total potentials 
of the reactants equals that of the products. The electrochemical electrode 
equilibrium is 


Ag + Cl'^ AgCl + e" (Al-4) 

where e“ stands for one mole of electrons in the silver metal. Thus the 
equilibrium conditions for the two electrodes respectively are 


★ 

The total potential of component i, is the sum of the contributions of 
the pressure-dependent, concentration-dependent, and electric field-dependent 
parts of the (chemical) potential (12, p. 94, 469, 480) 


yV + (p - 0.1013) v. + RTlna i + z.3^ 


(Al-1) 


where y9 is the standard potential, z, the valency of the particle (+1 for 
+ ‘ - * 

Ag , -1 for Cl , zero for species carrying no appreciable net charge, e.g. 

Ag and AgCl ) and <P the electrical potential (volt) at the locations considered: 


AgCl 


y AgCl + ( p " 0J013 ) v AgCl “ y AgCl 


(Al-2) 


and similarly for silver: 

5ftg = % i(AP3r 

Note that the pressure is measured in dekabar (0.1013 dekabar = 1 atm). 
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(Al-5) 


ij 11 + n 11 = ii** + n** 

y Ag y Cl“ y AgCl V 


y Ag + y CT y AgCl + V 


(Al-6) 


Subtract (Al-6) from (Al-5) 

Av M * AB C1- = 4M AgCl + 4i V (A, - 7 > 

V 

Because of the small volume of the (non-hydrated) electrons , the . 
effect of pressure on the total potential of the electrons is assumed to 
be negligible. Also, since the electron concentration in the wire, and 
hence y e _,are nearly uniform, it follows from (Al-1) that 

Ay e _ (5 Au e _ - 3fA4>_) - -TA4>. (Al-8) 

where A<f)_ is the electrical potential difference between the electrodes. 

Substitute for Ay e _ in eq. (Al-8) from (Al-7) 

A V = -^A<f»_ = Ay cl _ - (Ay AgC1 - Ay ftg ) (Al-9) 

According to "Poynting's equation" (12, p. 304), the variation with 
pressure of the free energy of one mole, G_./n.. of a pure component i, 
is related to the total volume, V.. , of mole of this component by: 

3G/n . \ 

1 = V./n. (Al-10) 

9p It 

Applying (Al-10) to pure AgCl and Ag respectively, we obtain: 

- - A ( G AgCl /n AgCl h = _ Ay AgCl = V AgCl _ Ap (A1 -1 1 ) 
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A ^ G Ag /n Ag^T 


(Al-12) 


= A % = v Ag A P 

Here we have used the fact that, by definition of the partial molar 
quantities, the molar free energies, 2^/n^ and the molar volumes V^/n^ 
of pure substances are identical with the partial molar free energies, 
vu , (chemical potentials) and partial molar volumes, v^ , respectively. 

Substituting Ay Ag ci and A y Ag from (Al-11) and (Al-12) respec- 
tively in (Al-9), we obtain 

A M e - = - - Ap(v AgC1 - V Ag ) (Al-13) 

Therefore the potential difference, A<j>_, between the two Ag/AgCl 
electrodes is 

A<f>_ = -(l/^)[Ay cl - - Ap(v AgC1 - V Ag )] = -(l/^)[Ay_ - Apv c] ] (Al-14) 

where the symbol v^ replaces (v AgC -| - v Ag ) and symbol y_ replaces 
U c ^-, for the sake of brevity. 

Comparison of this result with equation (12-6) p. 150 of ref. (2) shows 
that in the latter, the term Ap(v Ag Q-| - v Ag ) has been neglected. 

Let us examine to which extent this term affects the evaluation of 
the results of our experiments. To do so, we calculate first the numeri- 
cal magnitude of the expression Ap(v Ag ^-| - v Ag ): 

From ref. (13), p. 648, the specific gravity and molecular weight 

o 

of silver are 10.5 g cm and 107.88 respectively. Hence 

i: = 107 - 8 ? 3 10,28 - cm- -mole" 1 - . 

Ag 10.5 g cm"-* 
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and, from the same source 


- _ 143.34 on n „„3 T _-l 

v AgCl - — — " 25 

Hence - 

V C1 

= (25.8 - 10.28) cnf* mole"^ Ap 

0.965 x 10° amp sec eq - ^ ^ 

-15.5 


- 1 .61 X 10' 4 r— ~ Ap 
amp sec r 


For Ap = -0.2 dekabar (about 2 atm) this yields 


3? <v AgCI 


- V Ap = 1 


61 x 10 -4 x 0.2 


cm3 dekabar = 3.22 x 10- 5 volt 
amp sec 


i.e. -0.032 millivolt for 2 dekabar. 

The streaming potential for the C-103 membrane considered in this 
report was more than one order of magnitude higher. We conclude that as 
a very first approximation, eq. (12-6) of ref. (2) may be used, but that 
the correction term Ap(v Ag £.| - v Ag ) is not entirely negligible. 

Returning to eq. (Al-14), we substitute for Ay c i_ in terms of 
the definition of the total potential [eq. (Al-1)]: 


A<j,_ = -(I/SHCAuq-,- + v cl _Ap - ^A<j> - Ap(v AgC1 - v Ag )] 
= -(l/S^UPci- “3^ ” Ap ^AgCl " ^Ag " ^Cl - ^ 
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(Al-15) 


Now, if the volume change in the electrochemical reactions at the 
electrodes [eq. (Al-4)] is negligible, then the last term of equation 
(Al-15) is negligible (15b). 

Therefore : 

A4>_ - - C^viqi _/ 3^ ) + Acf> (Al-16) 

One can calculate the potential drop across the membrane, A<(>, 
which represents the difference in electrical potential between the 
two solutions in the immediate vicinity of the membrane surfaces. by 
simply adding a "concentration correction", (Ay^./^) = (RT/^)ln(aQ-|_/a^_) 
to the potential difference, A<J>_ , measured between the Ag/AgCl electrodes: 

A<f> * A<j>_ + (Ay^./Sf) (Al-17) 

A(J> appears as a component of the general driving force of the Onsager 
set of flux equations (3)-(5). Note that eq. (Al-17) was presented (without 
proof) in the 1968 Annual Progress Report. 


* , - - . 
Note that this term is likely to be much smaller than the term ApWa n - v^ ) 

discussed in connection with eq. (Al-14). 9 9 
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A. 2. Calculation of M-K Conductivity Coefficients, L 


All experimental values used are from Table 3. 


~ee 

R 

R 

RA 

P 


JA ? 2 


pd 

A 


1 


pd^ 

4 pd = RA 


0.28 ohm (from Figure 4) 

0.28 Q x 20.25 cm 2 = 5.6V ohm cm 2 
pd _ 5.67 


0.017 


= 332 ohm cm 


k = 1/p = 3.01 x 10~ 3 ohm - ^ cm~l 

(L eeWo = l/(pd^ 2 ) = 1/(5. 67 x 0.965 2 x 10 10 ) = 1.89 x 10~ 1] joule'Wsec'W 2 


L se ^ from eq * ( 42 ^ 


^ L se^Ay=0 ' t + L 


ee 


where 


t = (_§ ) 

c + " M‘/^ ; Ay=0 


^ L se^Ay=Q t+L ee 


-8.69 x 10" 9 mole cm" 2 sec"^ x 0.965 x 10^ coul eq"^ 

-0.971 x 10 3 amp cm -2 

= 1.64 x 10"^ eg 2 joule"^ sec~^ cm~ 2 


= 0.864 


'we 


From eq. (44) (L we ) Ay=0 = SAeeVo 
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where 


t w = cv* 1 ^ Vo 


•99.79 x 10 9 mole cm 2 sec" 1 x 0.965 x IQ 3 coul eg" 1 = g g2 

0.971 x 10~ 3 amp cm" 


^ L we^Ay=0 t w L ee 


9.92 x 1.89 x 10" 11 = 1.88 x IQ" 18 joule" 1 eg 2 sec" 1 cm" 2 


"ew 


From eq. (52) : 

< L ew>A U =o4 CyL ee (SPD + ¥ ! ' L «e 5 «3 


w 


From streaming-potential line,SPD = 0.176 x 10~ 3 r^~ x 9.87 atm 


SPD = 1.737 x 10" 3 volt dekabar 


atm 

-1 


dekabar* 


Ay=0 18.02 


1 moles cm -3 


0.965 x 10 3 coul mole" 1 


1 .898 x 10 


-3 


/ 3 -I** \ 

_ n o _ 7 j , 15.5 cm mole 

x 1.89 x 10 joule eq sec” cm. *1.737 x 10”' 5 volt dekabar" + g 965 x 10 5 coul mole" 1 


- 1.64 x 10 11 eq 2 joule" 1 sec~ 1 cm" 3 x 17.2 cm 3 mole" 1 ** 


( L ew) Ay=0 = TfW £ 348 x 10_11 " 28,2 x 10_11 ^ = 1-77 x 10~ 10 joule~ 1 mo1e 2 cm~ 2 sec~ 1 


1 dekabar = 10 bar = 9.8 atm. This unit was introduced (8 ) because it 
eliminates conversion factors when all electrical and osmotic measurements 
are expressed in units containing the practical wattsec = joule (1). 

From appendix A.l . 

icicie 

Ref. 21, p. 250, eq. (8-5-4). Parameters for equation from p. 253, Table 8-5-1. 


fl 


"ww 1 


ws 


= L , and L 
sw ss 


To calculate these coefficients (all at Ay i- 0), we first calculate 

c c 

Ay and Ay,,, 
s w 

c^ = l.OOlg x 10 4 mole cm 3 {= Q.lOOlgN) 
c^ = 0.495 4 x 10' 4 mole cm -3 ( = 0.0495 4 N) 

The mean activity coefficients, y +t of Na + (or Cl") are obtained by 

"Jc 

interpolation from ref. (22, Appendix 8.9, p. 466 ). Since the solutions are 
rather dilute, ratios of molar concentrations are sometimes substituted for 
ratios of molality throughout the following calculations: 


c' = 0.0495 x 10" 3 mole cm" 3 (m l = 0.0497 mole NaCl/1000 g H 2 0, vfif = 0.223); y+ = 0.821 
c" = 0.1001 x 10" 3 mole crrf 3 (m" = 0.1005 mole NaCl/1000 g H 2 0, M = 0.317); y+ = 0.780 

( a s /a s) = c 't+/(c"y;) = (0.0497 x 0.821 )/(0. 1006 x 0.780) = 0.52 

Ay^ = RTln(a;7a|) 2 = 2 x 8.13 joule mole'Vx )' 1 x 298(°K)[-ln(0. 52)] = 3169 joule mole~ 

0.65393 

The water activities are found by linear interpolation of published data 
(ref. 22, Appendix 8.3). 


Comparison with other' references showed -that -the- term. "activity _ coefficient of 
electrolyte" used in the heading of this table stands for the mean" activity ' 
coefficient, y ± , of the Na + ions (equal to the mean activity coefficient of the 
Cl" ions). The activity coefficient of the salt, Y^ a Q-j » equals y± 04* Table 1) 
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m = 0.0497 mole NaCl/1000 g H 9 0: a' = 0.998332 

L. W 

m = 0.1006 mole NaCl/1000 g H o 0: a" = 0.996626 

l w 


K = RTin < a w /a ;> ■ rt '"(4v ; w> 

= -8.13 joule mole -1 (deg K)" 1 x 298 °K x (0.001706/0.997) 

C . - *1 ★ 

Au = -4.15 joule mole 
w 


We wish to substitute numerical values in eq. (58). From 
Table 3, the dialysis «-* osmosis experiment 145, yielded the following flux- 
force ratio: 


/ J W 


\ -Ay, 

\ 


15.3 x 10 9 mole cm -2 sec 1 


w i=0 
Ap=0 


4.15 joule mole 


-1 


-9 -1 2 -2 -1 

3.7 x 10 joule mole cm sec 


* c c 

It is of interest to compare the ratio Au /Au to the corresponding value 

* ” c c 

calculated from the Gibbs-Duhem equation [eq. (9)]. We obtained (Au s /Au w ) = 
31 69/ (-4.15) = -764, whereas the Gibbs-Duhem equation yields -c $ /c w = 
-55.5/0.075 = -740. The discrepancy is due to the fact that, strictly, the 
Gibbs-Duhem equation is applicable only to very small concentration 
differences, dc. When used for finite concentration differences, c^ - c^, 
proper integration has to be performed, instead of using simply the average 
-- salt- concentration, c $ . 
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✓ 


r> 


Substituting this experimental value in eq. (58) and taking L 


es 


as equal to L sg , we obtain 


3,7 x 10-9, L . 7.; 7 X J0 ^ ° / 

^ 1.89 x 10" 11 


1.88 x 10' 10 x 1.64 x 10' 11 


"ws 


1.89 x 10 


-n 




1.76 x 10 


-9 


1 .627 x 10 


-10 


3170 

(-4.15) 

> 

-764 


(3.7 + 1.76) x 10“ 9 - 124.6 x 10" 9 = L, ni - 764 L 

ww 


V_- 


..... I 


WS 


-119.1 x 10 


-9 


Another relation between and L ws is obtained from eq. (59). From 
Table 3, we obtain the following hydraulic permeability values: 


/ J 


\ 


v | 

i=:.i 

and 


-6 -I 

y — ^ cm s - ec — x 9.87 atm dekabar -1 = 3.69 x 10 -6 cm dekabar -1 sec -1 

1 atm 


0.304 

0.8 


x 9.87 


3.75 x 1Q~ 6 cm dekabar -1 sec -1 

“6 i 

Average 3.72 x 10 cm dekabar sec 


The hyperfiltration effect was too small to be measured. Hence we take 
the composition of the permeate equal to that of the solution 


J s /J w = c s /c w = 10~ 4 mole cm" 3/55. 5 x 10"^ mole cm"^ = 1.8 x 10"3. 


J v - J w 5 w + J s ; s “ V 5 w + ’•« * ,0 ' 3 ^ 
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y 


Take v w = 18.02 cm 8 mole - ^ [calculated from eq. (13.60), ref. (23)] 
v s = 17.49 cm 8 mole"^ [calculated from eq. (13.54), ref. (23)] 

Then J = 18.05 J 
v w 

and 


\ 


w 


-Ap 


i=0 

Ay=0 



3 IT 
18.05 cm mole 


x 3.72 x 10 ^ cm dekabar - ^ sec 


- A ”l V 


= 2.06 x 10 _/ mole dekabar 


- 1 cm’ 2 


sec 


-1 


Substitute this value in eq. (59), simplified because Ap w = 0 = Ay g 

2.06 x 10' 7 = 17.2 L + 18.02 L - 1j88 x 10 (1.64 x 10" 11 x 17.2 + 1.77 x 10' 10 x 

WS m 1 .90 x 10" ' 1 \ 

3.48 x 10' 9 

2.06 x 10" 7 + 3.46 x 10" 8 = 18.02 L; + 1 7.49 L lie 

WW WS 

V ^ / 

2.41 x 10" 7 

Rewriting the two simultaneous equations for and L ws> we have: 

L , - 764 L = -1.191 x 10‘ 7 

WW WS 

L + 0.971 L = 1.33 x 10“ 8 

ww ws 

765 L ws = 1.324 x 10~ 7 -*■ L WS ( = L SW ) = ^* 7 ^ x joule~^mole 2 cnf 2 sec 

L... = 1.33 x 10‘ 8 - 0.017 x 10" 8 =1.31 x 10“ 8 joule" 1 mole 2 C ra" 2 sec" 1 ' - 

wW 


18.02) 

y 
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✓ 


r 


2.26 


L ss is calculated from eq. (56) as follows: 

From Table 3 we obtain for the salt permeability coefficient 


-0.0713 x 10 9 mole cm" 2 sec"^ 


-Avi 


i=0 

Ap=0 


-3170 joule mole 


-1 


2.25 x 10’^ 4 joule" ^mole 2 cm” 2 sec~^ 


Substitute this in eq. (56): 


( 10 


-14 


(1.64 x 10 


"ss 


1.89 x 10 




+ 1.73 x 10 


-10 (1 .64 x 10 ]1 ) x 1.77 x 10 10 /■ 

1.89 x 10' 11 


1.53 x 10 


-10 


2.26 x 10' 14 = L ss - 1.423 x 10" 11 + 2.54 x 10" 14 


It is seen that the dominant term in this equation is (L L /L ) = 

OC u 5 CC 

1.415 x 10”^ . In other words, relatively large errors in the salt diffusion 

-14 

flux, J s , which determines the first term on the left, 2.26 x 10 , affect 

the value calculated for L gs relatively little; the salt flux is primarily 
determined by the presence of an electrical potential difference across the 
membrane, rather than the "diffusion force", -Ay^. 

From the last equation 

L $s = 1.423 x 10" 11 joule"^ mole 2 cm" 2 sec"^ 

Table 4 in the main body of this report is a summary of the M-K conduc- 
tivity- coefficients ,_ L, calculated in this appendix. 

*If it were possible to effectively short the electrodes during the 
dialysis *-*■ osmosis experiment (A<j>_ = 0), one could calculate L $s by 
substituting the salt flux measured under these conditions into eq. (35), 


4.15 

3169 
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arid thus eliminate this difficulty, because the resulting equation no 
longer contains the large term ( L L /L ) 

SG GS GG 

(J s V=0 = L ss ( ‘ ip s ) + L sw ( - Au w ) 

Ap”=0 

It should be noted, however, that in this case the electrodes would 
have to be placed very close to the membrane faces, since the results of 
this measurement are expected to vary substantially with the distance of 
the electrodes from the membrane. As discussed in the introductory part 
of the section "Method for calculation of the L-conductance coefficients 
for the M-K transport equations", all the L-coefficients are to reflect 
membrane properties rather than solution properties. Therefore the method 
of evaluation was developed by a reasoning involving a thought experiment 
in which the electrodes are brought infinitely close to the membrane sur- 
faces. It was shown that many transport measurements in which this condi- 
tion is not met may be substituted, however, without materially affecting 
the results. The exceptions to this rule are the measurement of the electri 
cal resistance (in which the resistance for the conditions prevailing in 
the thought experiment was readily found, however, by extrapolating the 
straight-line plot of resistance vs. electrode distance to the distance 
corresponding to the electrodes touching the membrane faces), and the 
measurement discussed in this footnote. Here, too, a suitable extrapola- 
tion method might prove applicable in the future; in this case the current 
through the electrode "short" has to be measured. 
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